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Abstract 

Using the current world's largest data sample of 'i/'(3770) decays, we present results of a search for 
the non-DD decay ^^(3770) — > KgK^. We find no signal, and obtain an upper limit of (T('i/'(3770) 
KgK^) < 0.07 pb at 90% confidence level (CL). Our result tests a theoretical prediction for 
the upper bound on B(^(3770) KgK^) based on a charmonia-mixing model. 



*Current address: Universitat Bonn, Nussallee 12, D-53115 Bonn 



2 



The '?/'(3770), the hghtest charmonium resonance above open-flavor-production thresh- 
old, decays dominantly into a pair of D mesons. Non-DZ) decays of '?/'(3770), which are 
OZI-suppressed yy, have received much attention as they were observed to constitute an 
unexpectedly large portion (?«20%}of the total hadronic decay rate in earlier measurements 
summarized in 0. Recent results 0,01 do not conflrm a discrepancy of this magnitude, but 
due to the experimental uncertainties, the measurements are not yet conclusive. Exclusive 
non-DZ) decays of the ^^(3770) are expected as of any other charmonium state and have, in 
fact, been observed j^la01- The discrepancy between the hadronic and D-pair production 
cross sections will be clarifled by identifying other components of ?/'(3770) non-DD. 

This article focuses on the reaction ?/^(3770) — > K^K^. This flnal state is particularly 
interesting also in the context of the perturbative QCD "12% rule" (which relates J/ip and 
'ip{1S) branching fractions 0]), with the help of the mixing scenario ^3]; in which the 
ip{2S) and -0(3770) are considered to be mixtures of the 2^5*1 and l^Di states of charmonium. 
Rosner propounds a possible explanation of the "p — tt puzzle" using the phenomenon 
of charmonia mixing. In view of this, it is interesting to explore the nature and degree of 
mixing between the ip{2S) and ?/'(3770), particularly in the case of those modes which show 
signiflcant deviation from the 12% rule. The pseudoscalar-pseudoscalar (PP) flnal state 
KgKl is in this category: an average of recent results [1, [l^ indicates the ratio B{4'{2S) — >■ 
KgKl)/B{J/il> KgKl), to be as high as (29.9 ±3.0)%, substantially enhanced compared 
to the prediction of 12%. In a scenario that attempts to explain this enhancement, Wang 
and collaborators ^ predict (1.2±0.7) x 10"^ < i3(V^(3770) ^ K^s^l) < (3.8±1.1) x 10"^. 
The search presented here has sufficient sensitivit y to test the upper bound of this prediction. 

We use data collected by the CLEO detector [1^ operating at the Cornell Electron Stor- 
age Ring The sample corresponds to an integrated luminosity of 281 pb~^ of e^e~ 
annihilations at a center-of-mass energy ^/s =3773 MeV. The CLEO-c detector configu- 
ration features excellent efficiency and resolution for charged particles and photons within 
93% of the solid angle. The detector components critical to this analysis, rendering the dis- 
crimination of signal events amidst background, are discussed in the following. The charged 
particle tracking system consists of a low mass wire inner drift chamber (ZD) suitable for 
low track momenta, followed by an outer drift chamber (DR). These two devices measure 
charged track three-momenta with excellent accuracy and achieve a momentum resolution 
of ~0.6% at 1 GeV/c. The DR also measures energy loss that is used to identify charged 
tracks. The drift chamber is surrounded by a Ring Imaging Cherenkov Detector (RICH), 
followed by a Csl calorimeter (CC), where two regions in polar angle (measured with respect 
to the beam direction) are distinguished: barrel (|cos(6')| < 0.81) and endcap (|cos(^^)| > 
0.81). The CC allows the detection of photons with an energy resolution of 2.2% (5%) for 
photons with energy of 1 GeV (100 MeV) and, in combination with the tracking system, 
provides the basis for excellent electron identification. 

The event reconstruction and final state selection criteria proceed along the lines of 
CLEO's i'i^S) — > KgKl analysis [l2|- Our strategy is to reconstruct only a single Kg 
in each event and demand nothing other than an accompanying based on the following 
criteria. We reconstruct the K'^ using its decay to two charged pions, and thus require the 
events to have exactly two charged tracks. We have taken into account the effect of the 
small (~ 4 mrad) crossing angle between the e"*" and e~ beams by performing a Lorentz 
transformation of all the laboratory quantities to the center-of-mass frame. We impose 
standard track selection criteria based on the number of drift chamber hits and geometric 
acceptance. We use both charged particle ionization loss in the drift chamber (dE/dx) and 
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RICH information to identify the two tracks as pions which arc used to reconstruct the 
mesons. We define the parameters VYDuj = Li — Lj and PID2ij = erf — Cj, where 
Lij{i,j = p,K,7i with i ^ j) are the hkehhoods given by the measured Cherenkov angles 
of photons in the RICH detector compared with the predicted Cherenkov angles for the 
^th^ jth pap^icig hypothesis, and crjj(i, j = p, K, tt, e with i ^ j) are the ratios of the difference 
between the measured dE/dx and the predicted dE/dx values normalized to their standard 
deviations for the hypothesis. To use the RICH, we first require that the track 

momentum be above the RICH threshold of 0.6 GeV/c for the pion hypothesis. We then 
discriminate tt from p, K in the following manner: (a) If the RICH information is available, 
we require (PIDi;rp+PID27rp) < and (PIDi^i^ + PID2,ri^) < 0. (b) If the RICH information 
is not available, we require PID27rp < and PlD27rA' < 0. We further reject background 
from Bhabha (e+e~ — >■ e+e~) and two-photon {e^e~ — > 77) events by an electron veto of 
tracks that satisfy the condition for the ratio of the CC determined energy Eqc and the 
track determined momentum p as 0.92 < Eqq/p < 1.05 and have |(Te| < 3. 

The pair of charged pion candidates are kinematically constrained to come from a common 
vertex. We require that the reconstructed invariant mass of the two pions be within 10 MeV 
(Ri 3.2 standard deviations) of the nominal Kg peak. We reject background from non-X^ 
sources by requiring the measured flight path of the Kg candidate before its decay to be 
greater than 5 standard deviations {k, 5 mm) with respect to the interaction point. In 
addition, we require that the K^ candidates originate from the e+e^ interaction point by 
demanding that their distance of closest approach is within 5 standard deviations with 
respect to the interaction point. 

Frequently, K^ mesons will produce a shower in the CC. However, as this shower docs 
not have a measured energy corresponding to the energy of the we do not attempt to 
reconstruct it, but merely allow for its existence. In order to reject contamination from 
anticipated background events with energy from neutral particles other than the X^, we 
impose selection conditions as follows. We first find the K^ direction by looking opposite to 
the K^ direction, after having constrained the K^ to its nominal mass and also considered the 
effect of the finite crossing angle. We then require the energy of the shower associated with 
neutrals closest to the K^ direction to be less than 1.5 GeV. This cut further rejects QED 
background events of the type e^e~ — > 77 (where one of the 7 undergoes pair production). 
We also require the sum of the energy associated with neutrals outside a region around the 
direction of the defined by the angle between the position vectors of the K^ and the 
shower in consideration (^angie > 0.35 radians), to be less than 300 MeV. We do not include 
the energy from showers which have an energy below 50 MeV, which are frequently due 
to electronic noise, in this summation. This cut is very efficient in eliminating hadronic 
background events from the following sources: 

1. ir*°(892)i?° + charge conjugate (c.c.) produced in: 

(a) e+e" 7 ^(25), followed by the subsequent decay ^(2^) K*^{892)K^ + c.c, 

(b) e+e~^ K*^{892)K^ + c.c. via e+e" 7* (continuum). 

2. K^K^ produced in: 

(a) e+e~— > 7 1/^(23), followed by the subsequent radiative decay ip{2S) — > 7Xco,2, 

Xc0,2 -^5-^5' 

(b) e+e-^ V'(3770) ^ 7 Xco,2, Xco,2 ^ K^s^l 
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(c) e~^e — > KgKg via e~^e 7* (continuum). 

We define a "good" shower as one that has an energy profile consistent with being a photon 
and possessing an energy above 100 MeV, and require events to have no good showers 
associated with neutrals outside the K^^ region, and at most one good shower inside this 
region. 

Most events that contain one or more 7r° decays will be eliminated by the above cuts. 
However, for even better rejection, we find it useful to introduce an explicit 7r° veto with 
a lower photon energy requirement. To identify a 7r°, we require a pair of showers not 
associated with charged tracks to have their energy distribution consistent with a photon 
even if they overlap with nearby clusters. In addition, we require each photon shower 
candidate to possess at least 30 MeV (50 MeV) of energy for a barrel (endcap) photon and 
kinematically constrain the pair to the known 7r° mass. We further require the difference 
between the unconstrained and fitted 7r° mass, normalized by its resolution, to be < 3 and 
the 7r° momentum to be > 100 MeV/c. We reject all events that have any 7r° candidates 
meeting the above criteria. This is based upon our understanding of the basic topology 
of the event and our anticipated background interactions which produce ir^ mesons. Many 
background events will fail more than one of these cuts, but for best background rejection, 
all are necessary. 

The above selection conditions were optimized by studying a Monte Carlo simulation of 
events using the evtgen generator and a GEANT-based jT^ detector modeling program. 
We simulated events with a sin^ 9 angular distribution, where 9 is the angle between the Kg 
and the positron beam in the center-of-mass system, as is expected for a vector resonance 
decaying to two pseudoscalar mesons. We also included initial state radiation effects. 

We define a scaled energy variable for each event as the ratio of the Kg energy to the 
beam energy as i?^o/i?beam- Since the Kg mesons from the signal will be mono-energetic, we 
expect the Ej^o/ E]^eiim distribution to peak at unity for the signal. Based upon simulations 
(as seen in Figure [T]), we determine the signal region to be 0.98 < Ej^o/ E]^ea.m < 1-02. 

The relevant backgrounds studied are the QED sources (e+e^ — > 77, l^l^) and the 
hadronic sources: (a) K*^{892)K^ + c.c. where the K*^{892) decays into and K^, giving 
two neutral kaons in the final state, which could become a Kg and a K^, and (b) KgKg 
where one of the Kg tt+tt^, while the other Kg —>■ tt^tt^ and thus can mimic the signal. 
Two additional sources of background are identified and studied. One of them is DD 
production. Our Monte Carlo simulation study of a generic DD sample more than twice as 
large as the data indicates that this background is completely eliminated by our selection 
criteria. The second source of hadronic backgrounds taken into account is the more pervasive 
K'^Kl events originating from 4'{2S) produced at 3773 MeV in either the tail of or the 
radiative return to the ip{2S), the latter of which is the bigger contribution. These events 
peak at i^^^o/ii^beam = 0.977 as seen in Figure [H which uses a simulated sample of such 
events analyzed in an identical manner to the data. Unfortunately, these events cannot be 
eliminated fully in this analysis by using the total four-momentum constraint typically used 
in '?/'(3770) studies, as in this analysis we do not reconstruct the complete event. 

Figure 121 shows the i^^^o /i^beam distribution for events in data after the application of all 
of the above selection cuts. We observe 8 events inside the signal region, and an asymmetric 
background in the low sideband EksI E\,ea.m < 0.98. Simulation studies indicate that the 
events in the low sideband are not from backgrounds to the K^ sidebands but rather from 
multiple hadronic sources, dominated by contamination from 7 'iIj{2S)—>- 7 KgK^^, which may 
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FIG. 1: The scaled Kg energy (i?/^o/£'boam) distribution for the -0(3770) simulation samples cor- 
responding to the signal (solid histogram) and radiative return (dotted histogram) processes using 
arbitrary normalization. The arrows mark the signal region. 



well extend inside the signal region. In order to account for this, we estimate the number of 
events expected inside our signal region originating from this radiative return background. 
Our estimate of the amount of this contamination, obtained using the radiative return 
hadronic cross section at = 3773 MeV ^S] and i3(^(2^) ^ K^s^l) is 9.5±1.6 events 
inside the signal region. The uncertainty on the background estimate is found to be 17.1% 
by taking the quadrature sum of the relevant sources which include the uncertainties on the 
radiative return cross section Js^, the '?/'(3770) luminosity, the detection efficiency as obtained 
from simulations _(th e first six components listed in Table H)), B^ipi'^S) KgK^) [l^ and 
B{K'g — s> TT+TT^) |lj]. From such estimates, we find that the radiative return is a major part 
of the background, and it saturates our data in the signal region. In this analysis, we do 
not consider continuum subtraction from the resonance yield as this final state is forbidden 
to be produced via electromagnetic interactions under the SU(3) symmetry of fiavor fisl], as 
confirmed in 1^ 1. 



The detection efficiency for the signal channel is calculated to be 45.8% from Monte 
Carlo simulations. This efficiency includes a correction due to the difference in detection 
efficiencies determined between data and Monte Carlo simulation. The study uses e+e~ 
— > 70, (f) —>■ K^K^ events from the continuum channel in the same data sample as used 
in the analysis of '0(3770) — > K^K^. A possible momentum dependence arising from a 
difference in the momentum spectra of the mesons from the continuum and resonance 
channels is accounted for in the form of a systematic error component in the overall error 
attributed to the selection of 3.7%. A hardware trigger requiring that two tracks be found 
within 20 cm of the event vertex eliminated very long-lived mesons. The efficiency of this 
trigger for signal events was 73.6%, which is included in the calculated detection efficiency. 

Based on our observed signal yield and estimated radiative return background, we find 
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FIG. 2: The scaled energy (-E^o /-Ebeam) distribution for the ^(3770) K^K^ decay in 
281 pb^^ of data. The arrows mark the signal region. 



no evidence for a signal and calculate an upper limit on the number of signal events us- 
ing the Feldman- Cousins approach for an observed yield of 8 and a background esti- 
mate of 9 events which translates into an upper limit (90% CL) on the cross section of 
a{ip{3770) —>■ KgK^) < 0.06 pb. In order to include systematic uncertainties, we first con- 
sider the systematic uncertainty on the estimated background. We take the la lower value 
of this estimate and calculate the 90% CL upper limit on the number of signal events using 
the Feldman- Cousins approach. We subsequently include an overall systematic uncertainty 
of 5.5% from various sources listed in Table HJ and measure (t('?/'(3770) K^Kl) < 0.07 pb 
at 90% CL. 

In conclusion, using data collected by the CLEO detector at the ■?/'(3770) resonance, we 
have studied the exclusive non-DD '?/'(3770) decay to the PP two-body mode K^Ki^. We find 
no evidence of a signal in 281 pb~^ of data and report an upper limit on the production cross 
section, (j('?/'(3770) K^Kl) < 0.07 pb at 90% CL. This indicates that this mode is far 
from being able to account for any possible significant non- DD rate from the ^(3770) |3,||. 
In order to compare this result with the model prediction we translate our upper limit 
on the cross section into an upper limit on the branching fraction using the resonance cross 
section of the ^'(3770) reported by CLEO [J. We obtain i3(W3770) ^ K^s^l) < 1.17 x 10"^ 
at 90% CL, which improves upon the current upper limit [20| of < 2.1 x 10~^ (90% CL) by 
an order of magnitude. Our measurement is below the upper bound of the model prediction 
presented in [lO|. This may provide new insight into the nature of S- and D-wave mixing 
in charmonia and thereby help clarify the perturbative 12% rule. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the A. P. Sloan Foundation, 
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Source 


Systematic error (%) 


Trigger 


2.0 


Simulation statistics 


0.6 


selection 


3.7 


Tracking 


1.4 


Particle ID 


0.6 


finding 


3.0 


Luminosity 


1.0 


BiK'^S TT+TT-) [17] 


0.1 


Total 


5.5 



TABLE I: Summary of systematic errors. 



the National Science Foundation, and the U.S. Department of Energy. 



[1] S. Okubo, Phys. Lett. 5, 165(1963); G. Zweig, CERN TH 401, CERN TH 412 (1964), unpub- 
lished; J. lizuka. Progress in Theoretical Physics Supplement 37-38, 21(1966). 

[2] J. L. Rosner, Annals Phys. 319, 1 (2005), e-Print Archive: |hep-ph/Q411Q03| and references 
therein. 

[3] Q. He et al. (CLEO Collaboration), Phys. Rev. Lett. 95, 121801 (2005). 
[4] D. Besson et al. (CLEO Collaboration), Phys. Rev. Lett. 96, 092002 (2006). 
[5] N. E. Adam et al. (CLEO Collaboration), Phys. Rev. Lett. 96, 082004 (2006). 
[6] T. E. Coan et al. (CLEO Cohaboration), e-Print Archive: |hep-ex /0509030 (submitted to 
PRL). 

[7] G. S. Adams et al. (CLEO Collaboration), Phys. Rev. D 73, 12002 (2006). 
[8] W. S. Hou and A. Soni, Phys. Rev. Lett. 50, 569 (1983); M. E. B. Franklin et al. (Mark II 
Collaboration), Phys. Rev. Lett. 51, 963 (1983); W. S. Hou, Phys. Rev. D 55, 6952 (1997); 
Y. F. Gu and X. H. Li, Phys. Rev. D 63, 114019 (2001). 
P. Wang, X. H. Mo, and C. Z. Yuan, Phys. Rev. D 70, 077505 (2004). 
P. Wang, C. Z. Yuan, and X.H. Mo, Phys. Rev. D 70, 114014 (2004). 
J. L. Rosner, Phys. Rev. D 64, 094002 (2001). 

S. Dobbs et al. (CLEO Collaboration), e-Print Archive: |hep-ex/0603020| (submitted to PRL). 



[9] 
[10] 

[11] 

[12] _ 

[13] R. A. Briere et al. (CESR-c and CLEO-c Taskforces, CLEO-c Cohaboration), Corneh Univer- 
sity, LEPP Report No. CLNS 01/1742 (2001), unpublished. 
[14] G. Viehhauser et a/., Nucl. lustrum. Meth. A 462, 146 (2001); D. Peterson et al, Nucl. 
lustrum. Meth. Phys. Res. A 478, 142 (2002); M. Artuso et al, Nucl. lustrum. Meth. A 554, 
147 (2005). 

[15] D. J. Lange, Nucl. lustrum. Meth. Phys. Res., Sect. A 462, 152 (2001). 

[16] R. Brun et al, Geant 3.21, CERN Program Library Long Writeup W5013 (1993), unpubhshed. 
[17] S. Eidelman et al (Particle Data Group), Phys. Lett. B 592, 1 (2004). 



8 



[18] H. E. Habcr and J. Pcrricr, Phys. Rev. D 32, 2961 (1985). 

[19] G. J. Feldman, and R. D. Cousins, Phys. Rev. D 57, 3873 (1998). 

[20] M. Ablikim et al. (BES Collaboration), Phys. Rev. D 70, 077101 (2004). 



9 



